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Abstract 

The corrections of two-photon exchange on deuteron electromagnetic form factors are estimated 
based on an effective Lagrangian approach. Numerical results for the form factors Gc,ai,q of the 
deuteron with the corrections are compared to its empirical ones. Moreover, the two new form 
factors, due to the two-photon exchange, are analyzed. Possible way to test the two-photon exchange 
corrections to the deuteron form factors is discussed. 
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1 Introduction 



We know that the electromagnetic (EM) form factors of the proton and deuteron are usually extracted 
from the measurements of the differential cross sections of ep and eD elastic scatterings and from the 
Rosenbluth separation method [1], which is based on one-photon-exchange (OPE) approximation. For a 
long time, the extracted Q^-dependences of the nucleon EM form factors are believed to behave like a 
simple dipole form. For the proton electric and magnetic form factors, j^^, one conventionally assumes 

G^(Q2) = G^(Q2)/^^ ^ 1/(1 + Q\GeV^)lQ.nf, (1) 

where ^.p — 2.79 is the proton magneton. Recently, the new experiments of the polarized ep elastic scatter- 
ing were precisely carried out at Jefferson Laboratory [5] . The polarization transfer scattering experiments 
oie+p e+pshow that the ratio RP = ^ipGl.{Q^)/Glj{Q^) behaves like RP{Q^) - 1 - 0.158Q^. It 
means that is no longer a simple constant as implied in eq. (1). It monotonously decreases with the 
increasing of Q^. 

One way to resolve this discrepancy, at least partially, is to take the effect of the two-photon-exchange 
(TPE) into account [3 H E [S] . Usually, it is believed that TPE is strongly suppressed by EM 
coupling constant asM (~ 1/137). However, it was argued [5] that due to a very steep decreasing of 
the nucleon EM form factors, the TPE process, where the is equally shared by the two exchanging 
photons, may be compatible to the OPE one. Some calculations of the TPE corrections to the ep elastic 
scattering have been done recently [H HI [3 [HI [3 [5] , where only the nucleon state is considered as an 
intermediate state. The calculations were extended further with other nucleon resonances, like A, Pn 
and Di3 states, being considered as the intermediate states [10]. There were also several other works 
about the TPE effect on the proton charge radius and on the parity-violating [HI [T2] in the ep scattering. 
The effect on the EM form factors of the nucleon in the time- like region was estimated in Refs. [T^ [Tl] . 
According to the analyses for the TPE effect on the nucleon EM form factors in the literature, it is known 
that the TPE corrections not only modify the conventional nucleon electric and magnetic form factors, 
but also provide a new form factor, 12^, to the nucleon. 

The TPE corrections to the deuteron (spin 1 particle) EM form factors and to the e+ ~* D + D 
process have been also discussed in Refs. [TSJ HH HZ] qualitatively. In analogy to the TPE effect on the 
proton EM form factors, TPE not only modifies the conventional three EM form factors of the deuteron, 
but also provides new form factors with new structures. The general discussion of the structures of the 
three new form factors can be seen [T5J [TO] . We know that the deuteron is usually regarded as a weekly 
bound system of a proton and a neutron (see Fig. 1). Many calculations for the EM form factors of the 
deuteron, with the OPE approximation, have been performed in different approaches in the literature 
(see for example, Refs. [H] [H] [201 [2I])- Recent calculations based on an effective Lagrangian approach 
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Fi gure 1; Dcuteron mass operator. 

[551 [131 have shown that this approach can reasonably explain the deuteron EM form factors with phe- 
nomenological including two-body operators. 

To study the TPE effect on the deuteron system in our effective Lagrangian approach, we note that 
the deuteron EM form factors receive the TPE corrections from three different sources. The first one (see 
Fig. 2) is that the two photons directly couple to the contact points (the contact point A (or B) of Fig. 
f is the one connects the deuteron to its composites) . The second is that one of the two photons directly 
couples to one of the nucleons and another to one of the contact points (see Fig. 3). The last one is 
that the two photons respectively couple to the two nucleons (see Fig. 4). It has been proved that gauge 
invariance preserves in our effective Lagrangian approach only when the three kinds of the two-photon 
exchange diagrams are considered simultaneously [24j . 

In our previous work [51] , only a part of the third type of the TPE corrections to the EM form factors 
of dcuteron is considered, where the TPE corrections to the EM form factors Ge.m and to I27 of the 
proton and neutron are directly employed to study the deuteron properties following the formalism of Ref. 
[9]. In the approach, only one new form factor appears. The TPE effect considered in [25] is represented 
by Figs. 4(a), 4(b) and their cross-box diagrams. In this paper, to extend the work of [^ further, we'll 
simultaneously study the three sources of the TPE effect in Figs. 2-4. It should be stressed that although 
the contribution of the coupling of a photon to the contact point is expected to be smaller than the one of 
direct couplings of the photon to the nucleons, this type of couplings is needed in order to guarantee gauge 
invariance. This paper is organized as follows. In section 2 the above mentioned two-photon-exchange 
effect in the eD elastic scattering is briefly discussed. Numerical results and conclusions arc given in 
section 3. 
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Figure 2: Diagrams for the first type of the two-photon exchange effect. The cross-box diagrams are imphed. 

2 Two-Photon-Exchange in the eD elastic scattering 

According to the OPE approximation, the electromagnetic form factors of the deuteron are defined by 
the matrix element of the electromagnetic current J^{x) 



< p'c, A' I J,,(0) I PD, A >= -eoi Gi(Q^)r (A') • e(A) 



(2) 



-G3(g 



2,(e'*(A')-9)(e(A)-g) 



p, + G^iQ') um'*{^') ■ q) ~ C(A')(e(A) ■ q) 



where p'^,, A' (or po, A) denote the momentum, helicity, and polarization vector of the final (or initial) 
deuteron, respectively. In eq. (2) q = p'^ — pjj is the photon momentum, P = pjj + p'^, = —q^ is the 
four-momentum transfer squared, Md is the deuteron mass, and eo is the charge of the deuteron. In the 
one-photon exchange approximation or Born approximation, the unpolarized differential cross section of 
the eD elastic scattering, e(fci, si) -I- D{po,S,) — » e.{k'i, S3) -f- D{p'jj,£^'), in the laboratory frame is [55] 



da 
dJ7 



da 

dh 



loiOPE), 

■t 



(3) 



where 9 is the scattering angle of the electron, {da/dfl)Mott is the Mott cross section for a structure-less 
particle with recoil effect, and the two structure functions are 

2 



MQ') = GUQ') + -TDGi,iQ') + -Tt,Gt^iQ'), 
B{Q') = T^rD{l+rD)Glj{Q^). 



(4) 



In eq. (4) tu = (5^/4M^, and Gi\i, Gc and Gq are the deuteron magnetic, charge and quadrupole form 
factors, respectively. They can be expressed, in terms of Gi, G2 and G3, as 



Gm — G2, Gq — Gi — G2 + (1 + T£))G3, Gc — Gi + —tdGq. 



(5) 
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Figure 3: Diagrams for the second type of the two-photon exchange effect. The cross-box diagrams are implied. 

The normalizations of the three form factors are Gc(0) — 1, Gi\i{0) = 1.714, and Gg(0) = M^^Qo = 
25.83. Note that in eqs. (3) and (4), there are two unpolarized structure functions A and B, and three 
independent form factors Gc, Gq and Gm for the deuteron. To determine the three form factors com- 
pletely, one needs, at least, one polarization observable. The optimal choice is the polarization T20 (or 



Considering both OPE (C = —1) and TPE (C ~ +1), and taking Lorentz, party, and charge- 
conjugation invariance into account, one obtains the most general form of the eD elastic scattering 

2 6 

•^''^ = S3)7M^(fci, Si) E ^^^^^ (6) 

^ i=i 

where 



and 



11 _ 



Mf = 



iC ■ OP"", 

enr • 9) - (e • 9K'*'^' 



(7) 
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Figure 4: Diagrams for the third type of the two-photon exchange effect. The cross-box diagrams are implied. 



■ K) + ■ K)C*'' 

1 



(8) 



where K = ki + k[. General speaking, the form factors G'^, with z = 1,6, are complex functions of 
s = (pd + and = — (fci — k[Y . They can be expressed as 



G[[s,Q^)^G,{Q^) + Gf\s,Q^), 



(9) 



(2) 

where Gi corresponds to the contributions arising from the one-photon exchange and G\ stands for the 
rest which would come mostly from TPE. In the OPE approximation, G'^ — G'^ ^ G'q = 0. It is easy to 
see that Gi {i = 1, 2, 3) are of order of {aEAi)^ and Gp' {i = 1, ...6) are of order uem- 



To consider that a deuteron is a weakly bound state of a proton and a neutron, we take the following 
effective interaction between the deuteron and its composites (pn) }23j 

Cd = goD^+ix) J dy'fDiy^)p{x + ^y)Gj^n{x - ^y) + H.c, (10) 

where C is the charge conjugate matrix, D^, p and n are the fields of the deuteron, proton and neutron. 
The correlation ftmction $£i in eq. (10) characterizes the finite size of the deuteron as a pn bound state 
and depends on the relative Jacobi coordinate y, in addition, x being the center-of-mass (CM) coordinate. 
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The Fourier transformation of the correlation function reads 

^D{y^) = l^^~'''^^i~p')- (11) 

A basic requirement for the choice of an exphcit form of the correlation function is that it vanishes suffi- 
ciently fast in the ultraviolet region of Euclidean space to render the Feynman diagrams ultraviolet finite. 
Here, we adopt a Gaussian form ^^{p^) = exp{~p'^/ A'j^) for the vertex function, where pE is the Eu- 
clidean Jacobi momentum of the deuteron, and is a size parameter. It characterizes the distribution 
of the constituents inside the deuteron. 

We know that the low energy theorem [29] provides a model independent test for the reliability of 
different approaches [30j . For the photon deuteron (spin-1) Compton scattering, the low energy theorem 
has been discussed extensively in the past ^31J. A complete treatment on this issue has been given by Ref. 
[32j . To the first order of the photon energy lo, the low energy theorem tells that the forward Compton 
scattering amplitude off the deuteron target is [30] 

2 2 

^''^ " "7^'' ■ '~ ' ^^'^ • (e^' X 6) + 0{lu'), (12) 

with e* (or e '), S and fio being the initial (or final) photon polarization, the deuteron spin and its mag- 
netic moment in unit of e/2M£), respectively. In eq. (12) the first and second terms are the Thomson 
and the spin- flip ones. The latter is proportional to the deuteron anomalous magnetic moment squared 
k|) = (M-D — 2)^ and associates to the well-known Drell-Hearn-Gerasimov sum rule of the deuteron (see 
[30|, I33j for example). In our effective Lagrangian approach 23 , the effective current of photon-deuteron 
has the correct structures like eq. (2) and the numerical calculation shows that the obtained magnetic 
moment of the deuteron /ijj is around 1.7 (in unit of e/2AI]j), which reasonably agrees with the exper- 
iment data. Consequentially, it is expected that the forward Compton scattering amplitude based on 
our effective approach is consistent with the low energy theorem. It should be mentioned that the above 
correlation function of eq. (11), in the non-relativistic approximation, stands for the wave function with 
only S-wave of the deuteron, which does not contain any D-wave component. To reasonably explain 
the data for the deuteron quadrupole moment, we have to phenomenological include two-body operators 
[23j . A detailed comparison of the photon-deuteron Compton scattering amplitudes in the low photon en- 
ergy region of our approach and of the low energy theorem will be explicitly given in a our separate paper. 

In our approach, the coupling gu of < p^,A | pn >= gui *(A) is determined by the compositeness 
condition [24l [34] [35l [36l [37| . It implies that the renormalization constant of the deuteron wave function 
is set equal to zero: 

Zd = 1- ^'oi^lh) = 0. (13) 
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Here, 

(14) 

is the derivative of the transverse part of the mass operator S^'', which conventionally splits into the 
transverse T^d and longitudinal S|^, parts as: 

<^5fs,,b|,) + ^S^bi), (15) 
Pd 

where g'^^ — g"'^ ~ p'^p^ /p^ , and g'^Pa = . The mass operator of the deuteron in our approach is 
described by Fig. 1. If the size parameter is fixed, the coupling is fixed too according to the 
compositeness condition (13) (see detail in [23]). Here, we reiterate that since Figs. 2, 3 and 4 are taken 
into account simultaneously, gauge invariance preserves in our effective Lagrangian approach. 



3 Numerical results and conclusions 

To proceed a numerical calculation, we adopt the parametrization forms of the nucleon EM form factors 
given by Mergell, Meissner and Drechsel 38J. Here we follow the numerical technique of [39, to simplify 
one of our loop integrations. The loop momentum of the box- type Feynman amplitude is parametrized 
in a such way, that the denominators of Green function are (=fk + q/2)^ for the two photons (see Fig. 5 
where the cross-box diagram is explicitly shown) , whereas for the electron (e) and the constituent nucleon 
(N), they have the forms of (e) = {-n + fCY - ml and {p) = (±/C + VY - Mjj with 

IC = l{h+k[) = lK, V=\{p + p'). (16) 

The sign "-"("-|-") is for the direct (cross-box) diagram in Fig. 5. Here, it should be mentioned that 
the assumption of 39J means that each of the photons carries approximately half of the transferred 
momentum q. It is justified on the bases of Ref. [Sj. Moreover, the assumption also means that a rapid 
decreasing of the form factors is employed such one can neglect the dependence on the loop momentum 
K in the denominators of the photon Green function as well as in the arguments of the form factors. 
This results in ultraviolet divergences of the loop momentum integrals with respect to the momentum k. 
Thus, a step function 9{M'^T]\[~ | |) is introduced in the loop integration. It is equivalent to apply a 
cut-off restriction | |< Mf^r. For the photon Green function, we have 



with tm = TaW- calculation for the effect of TPE based on the effective Lagrangian of eq. (10), we 

face two loop-integrations. One is the loop integration with respect to the intermediate momentum of k. 
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n n 

(a) (6) 
Figure 5: Feynman diagrams for two-photon exchange: box diagram (a) and crossed box diagram (b). 

and another is the one with respect to the intermediate momentum k of the composites of the deuteron 
(see Fig. 5). To simpHfy the numerical calculation further, we also use the soft approximation for the 
integral variable k, in the first loop integration. 



Based on the above assumption and considering the TPE effect shown in Figs. 2-4, we may estimate 
the TPE corrections to the deuteron EM form factors in the present effective Lagrangian approach. The 
effective EM interaction Lagrangians have already been given explicitly in Refs. [531 121] • With those 
Lagrangians, we can correctly get the normalization conditions for Gc(0) and Ga/(0). In Figs. 6-10, we 
plot our numerical results for the contributions of the TPE effect to the deuteron electromagnetic form 
factors of Gc, Ga/, Gq and to the two additional form factors G5 and Gg. Two different scattering 
angles, 9 being 7r/2 and tt/IO, are selected in order to check the 9 dependences of the observables. In 
Figs. 6-8, the ratios stand for 

nc.M.Q = -7^. — 777«T' V-^) 

where g''^\j q{s , Q"^) represent for the TPE contributions. The individual contributions of Figs. 2, 3, 4 
and their sum to the deuteron form factors are shown explicitly. G^^ g in eq. (18) are estimated by 
the parametrizations of Ref. |40| (^-independent form) as the empirical data. The two maximum points 
in Figs. 6 and 7 are due to the two crossing points of the charge Gc and magnetic Gm form factors at 
about Qlrossing ~ 0-5 GeV^ and Qlrossing 2.0 GeV^. Here, different from the form factors of the 
nucleon, the form factors of the deuteron have the crossing points. From Figs. 6-8, one cannot explicitly 
see the 6'-dependence of the three ratios, since the dependences are strongly suppressed due to the fact 
that the denominators of the ratios in cq. (18) are 6'-independcnt. However, the ^-dependences can be 

(2) (2.) 

seen explicitly in Figs. 9-10 for the two new form factors Gg = Gg and Gg = Gg . It should be men- 



9 



tioncd that we do not have the extra form factors G4 contributed by the TPE effect as shown in eqs. (6) 
and (8) since we adopt the assumption of [39^ , where the K-dependence terms in the numerator are ignored. 

In our calculation, we have one parameter Aj^ in the correlation function. According to the condition 
that the deuteron is bound as <| r^^ |>< 0.02 GeV^ [T^, we select a typical value for the parameter: 

~ 0.30 GeV which is consistent with the one used in Refs. [531 US]- Our estimates for the ratios of 
the deuteron electromagnetic form factors of GJ^ q tell that the TPE effect is small. To analyze the 
contributions of Figs. 2, 3 and 4, one sees that in the low region, the contributions of Figs. 2 and 3 
are smaller than the one of Fig. 4. When increases, the contributions of Figs. 2, 3 and 4 increase too. 
Moreover, the contributions of Figs. 2 and 3 to the form factors of Gc.q are always smaller than that of 
Fig. 4, whereas the one of Fig. 2 to Gm becomes compatible to the contribution of Fig. 4. Since the total 
contributions of Figs. 2, and 3 to the three conventional EM form factors of Gc,m,q are very small, and 
the 6'-dependences of the ratios from the TPE corrections are suppressed, it is not easy to directly test 
the TPE effect from the three form factors. 



However, it is expected that one may test the TPE effect from the polarizations of deuteron, since the 
TPE effect is 6'-dependent and the obtained new form factors G^^ are 6'-dependent too. Consequently, it 
is reasonable to find the TPE effect in some polarizations and particularly in some angle limit. We know 
that if only one-photon-exchange is considered, the double and single polarization observables are 



P. 



P. 



-TD--;-i^(^n-GMGQ, 
\ ^ ^^d{td + l)tan^ ^ G\i . 



(19) 



Clearly, these two polarizations become vanishing when 9 is very small since they are and tan^^- 

dependent, respectively. However, when the TPE effect is considered in the small angle limit, its contri- 
bution is 



SP^. 



2TliC0t- 



^TD + 1 



Ga i?e(G^) 



Gi -4G2 + 2(Ti3 + l)G3)i?e(G[i) 



(20) 



and 



5P, 



2td 
3 V 



TD 



3 + 2{td + l)tan^ - ) G2 i?e (G^; ) 



'2{td + l)G2Re{G'^) 



(21) 



One sees that the TPE corrections to the polarizations do not vanish in the limit oi 9 ^ 0. In Fig. 
11, we display the ratios R{Pxz) = SPxz/Pxz for Pxz^ and R{Pz) — SPz/Pz for P^ calculated from eqs. 
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Figure 6: | Rc \ (%) for 8 = n/lO (a) and for 8 = n/2 (b). The dotted, dashed, dotted-dashed and solid curves represent 
the contributions from Figs. 2, 3, 4 and their sum. 

(19)-(21). The ratios should behave as l/ian^(|). One sees that the contributions from Figs. 2 and 3 
are found to be smaller than that of Fig. 4. Moreover, one finds the sizeable effect of the two new extra 
form factors, due to TPE, on the polarization observables Pxz and P^. The remarkable 6'-dcpcndcnces of 
the ratios are also displayed in Fig. 11. Therefore, a precise measurement of the deuteron polarizations 
in the small angle limit is expected to test the TPE effect. Since the deuteron form factors have crossing 
point Qlrossing^ ^^^^ cxpectcd to easily find the TPE effect at about Q"^ ^ Q crossing- 

Ref. [25] is the first one to numerically estimate part of the TPE corrections to the deuteron form 
factors based on our effective Lagrangian approach. In that work, the TPE corrections, to the EM form 
factors of the proton and neutron following the formalism of Ref. [9j , are simply employed to study the 
deuteron case. The corresponding TPE effect on the deuteron is shown by Figs. 4(a), 4(b) and their 
cross-box diagrams. Comparing the present results to those of Ref. [25 , one concludes that all the possi- 
ble TPE corrections are considered in this paper. Therefore, the present work gives a more systemically 
and sophisticated study of the TPE effect on the deuteron. Moreover, we directly calculate the TPE 
exchange effect with the assumption of [39 in this paper. Clearly, the present calculation gives more 
information about the new deuteron form factors since we predict form factors of Gg g simultaneously. 
The obtained results for the TPE effect are consistent with the ones of [25] qualitatively. Finally, one 
still cannot get any information about G'^, this is due to the approximate methods we employed here to 
simplify our numerical loop integration. 

To summarize, we are the first to estimate all the TPE corrections, as shown in Figs. 2-4, to 
the conventional form factors of the deuteron, Gc.m.q and of Ggg. Our numerical results of the TPE 
contributions tell that G(j q are small (less than 1%). However, G'^ q are clearly ^-dependent. The two 
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Figure 7: ] Rm \ (%) for 8 = vr/lO (a) and for = tt/2 (b). Notations as Fig. 6 



additional form factors are expected to be tested in the future measurements of the double and single 
polarization observables of Pxz{T2i) and Pz (Tin) hi the small angle limit and at about ^ Q crossing- 
Further work for an exactly full calculation of the two-photon exchange effect on the deuteron system, 
without using assumption of (39j . is in progress. 

Finally, this work is also designed to effectively treat direct electromagnetic interactions to quarks. 
It should be addressed that the present investigation of the two photon exchange mechanism recalls a 
new study of Compton scattering and it is shown that the local two-photon coupling to the same quark 
provides a fixed Regge singularity at J=0 ^T]. This subject is beyond the scope of the present work. 
However, it is of a great interest to see the issue for the deuteron target in our future work. 
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Figure 8: \ Rq \ (%) for 6» = tt/IO (a) and for 6 = 7r/2 (b). Notations as Fig. 6 
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